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Abstract 

Results are presented on if)(2S) and J /if) hadronic decays to KgpK^n and KgpK + n final states 
from data samples of 14 million ip(2S) and 58 million J/ip events accumulated at the BES II de- 
tector. No 0(1540) signal, the pentaquark candidate, is observed, and upper limits for B{ip{2S) — > 
69 -» K° s pK-n+K s pK + n) < 0.84 xl0~ 5 and B(J/%b ^69^ K° s pK-n+K° s pK + n) < 1.1 xlO" 5 
at the 90% confidence level are set. For single 0(1540) production, the upper limits determined 
by our analysis are also on the order of 10 -5 in both i^{2S) and J/ijj decays. 



PACS numbers: 13.25.Gv, 14.20.Jn 



I. INTRODUCTION 



Recently the LEPS Collaboration at Spring-8 reported the 4.6 a discovery of a new 
S = +1 state, the 0(1540), with a mass of 1.54 ± 0.01 GeV/c 2 and a width of less than 
25 MeV/c 2 , close to NK threshold in the reaction 7 12 C — > K + K~X [1]. Subsequently, the 
DIANA Collaboration at ITEP, CLAS at Jefferson Lab, and SAPHIR at ELSA claimed this 
narrow state, the candidate for a pentaquark state (uudds), in n + K + or p + K° decay 
configurations - all published in 2003 [2-4] . 

The work of LEPS was motivated, in part, by the work by Diakonov, Petrov and Polyakov 
who studied anti-decuplet baryons using the chiral soliton model Q]. In their model, the 
anti-decuplet was anchored to the mass of the Pn(1710) nucleon resonance, giving the 
pentaquark state + (spin 1/2, isospin 0, and S = +1) a mass of ~ 1530 MeV/c 2 and a 
total width of less than 15 MeV/c 2 . In particular, it is predicted to be an isoscalar. 

There are many other theoretical works to try to explain the properties of the 0(1540) 
with various quark models [6-9] or alternative approaches ^(|. Isospins and 1 are both 
possible; isospin 1 would lead to three charge states 9°, + , ++ . A reanalysis [11] of older 
experimental data on the f^ + -nucleon elastic scattering put a more stringent constraint on 
the width to be To < 1 MeV. Since the mass of the is larger than the sum of the masses 
of the nucleon and kaon, it is not easy to understand why its width should be so narrow, 
unless it has very special quantum numbers. Capstick, Page and Roberts point out that 
the + could be a member of an isospin quintet with charges from —1 to +3 where the 
Q = +3 state has a (uuuus) quark-model configuration. Decays of an isotensor + into 
pK° and pK + are isospin violating; hence an isotensor + is expected to be narrow jf|. 
The analyses of CLAS and SAPHIR support that the + is isoscalar, but the statistics of 
present experiments are limited. 

There are experimental questions concerning the 0(1540). For the four experiments [1- 
4], the exact shape of the background is very difficult to estimate since the + is close to 
NK threshold, and strong cuts have been applied to the event samples in all cases. The 
calculation of Dzierba et al. shows that kinematic reflections of meson resonances could 
well account for the enhancement observed in the K + n effective mass distribution at the 
mass of the purported + . They suggest that further experimental studies will be required 
with higher statistics, including varying the incident beam momentum and establishing 



the spin and parity, before claiming solid evidence for a S = +1 baryon resonance. It is 
important to understand the properties of the through systematic studies by different 
experiments. Compared with the above experiments, the data accumulated at the e + e~ 
collision experiment BES are relatively clean and have less background; therefore it is mean- 
ingful to investigate the pentaquark state G with the hadronic decays of the charmonium 
states ip(2S) and J/tf>. 

In this work, we search for the pentaquark state 0(1540) in ip(2S) and J /if) decays to 
K s pK~n and K s pK + n final states using samples of 14 million if>(2S) and 58 million J /if) 
events taken with the upgraded Beijing Spectrometer (BES II) located at the Beijing Electron 
Positron Collider (BEPC). These processes could contain decays to K s p, K + n (uudds) 
and decays to K s p, K~n {uudds). 

II. BES DETECTOR 

BES II is a large solid-angle magnetic spectrometer that is described in detail in Ref . . 



Charged particle momenta are determined with a resolution of a p /p = 1.78%y 1 + p 2 (GeV 2 ) 

in a 40-layer cylindrical drift chamber. Particle identification is accomplished by specific ion- 
ization (dE/dx) measurements in the drift chamber and time-of-flight (TOF) measurements 
in a barrel-like array of 48 scintillation counters. The dE/dx resolution is <JdE/dx — 8.0%; the 
TOF resolution is otof = 180 ps for Bhabha events. Outside of the time-of-flight counters 
is a 12-radiation-length barrel shower counter (BSC) comprised of gas proportional tubes 
interleaved with lead sheets. The BSC measures the energies of photons with a resolution of 



<Je/E ~ 21%/ J E(GeV). Outside the solenoidal coil, which provides a 0.4 Tesla magnetic 
field over the tracking volume, is an iron flux return that is instrumented with three double 
layers of counters that are used to identify muons. 

In this analysis, a GEANT3 based Monte Carlo simulation package (SIMBES) with de- 
tailed consideration of detector performance (such as dead electronic channels) is used. The 
consistency between data and Monte Carlo has been checked in many high purity physics 
channels, and the agreement is quite reasonable. 





III. EVENT SELECTION 



In ip(2S) or J/ip decays to K s pK~n and K s pK + n, the anti- neutron and neutron are 
not detected. The first step in the analysis is to select in four prong events the 7r + 7r~ pair, 
which composes the K s . Next the missing mass is calculated according to energy- momentum 
conservation. The events with missing mass close to the n(n)'s mass are selected. We use 
the same criteria and treatment for both ip{2S) and J/ip data. 

The first level of event selection requires two positively and two negatively charged tracks 
with less than ten neutral tracks. The charged particles are each required to lie within the 
acceptance of the detector and to have a good helix fit. 

The Kg meson in the event is identified through the decay Kg — > it + it~. For the selection 
of the Kg, both the 7r + and 7r~ are required to be consistent with being pions, namely, 
Prob P id(n) > 0.01. The definition of Prob P id is 

2 _ 2 1 2 

Xcom — XTOF ~t XdE/dx 

Prob ptd = Prob( X 2 com ,2), 

where Xtof an d XdE/dx are determined from the measured and expected TOF times and 
ionization information for the particle hypothesis of interest. If only TOF or dE/dx infor- 
mation is available, Prob P i d = Prob(x 2 , !)■ The four charged tracks can be grouped into a 
maximum of four possible 7r + 7r~ combinations. The combination with the invariant mass 
closest to the mass of K° s is chosen as the K s candidate, and the intersection point of the 
two tracks is regarded as the secondary vertex. Candidate events are required to satisfy 
\M n + n - — M K o \ < 15 MeV/c 2 , where M n + n - is calculated at the K s decay vertex. For the 
proton (or anti-proton) and kaon, we require Prob P i d {p) > 0.01 and Prob P i d (K) > 0.01, 
respectively. 

Candidate events are kinematically fitted (one constraint fit) under the assumption of 
a missing n(n) to obtain better mass resolution and to suppress the backgrounds. As an 
example, Fig. 1 shows the XxOpK-n distribution of i>{2S) — > K s pK~n candidate events 
after particle identification. We require x\o pK ^ n or x\v p K+n < ^ anc ^ use a further cut on 
the Kg decay length, L xy > 3 mm, to remove remaining backgrounds. 

Fig. 2 shows the n missing mass distributions of data and Monte Carlo for ip{2S) — > 
K s pK~n after the above cuts. The surviving events fall mostly within the range 0.9 - 1.0 
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FIG. 1: x 2 distribution of ip(2S) — > KgpK~n candidate events after requiring particle identifica- 
tion. The solid histogram is for data, and the dashed histogram for Monte Carlo simulation (phase 
space), where the Monte Carlo is normalized to the data in the first bin. 

GeV/c 2 , and we further require 0.9 GeV/c 2 < m missing < 1.0 GeV/c 2 . The missing mass 
distributions for the other modes if>(2S) — * KgpK + n and J /if) — > K^pK'n, K^pK + n are 
similar to those of Fig. 2. 

IV. ANALYSIS RESULTS 

After the above requirements, the individual mass distributions of ip(2S) — > K^pK'n 
and K$pK + n are shown in Fig. 3, and the scatter plot of K~n (Kgp) versus K$p (K + n) 
for ip{2S) — > K^pK~n + KgpK + n modes is shown in Fig. 4. No clear signal is observed 
in Fig. 3, which contains 19 (KgpK~n) and 10 (KgpK + n) events, or Fig. 4. We determine 
an upper limit for B(ip(2S) — > 00 — > K®pK~n + K^pK + n). The signal region is shown 
as a square centered at (1.540, 1.540) GeV/c 2 in Fig. 4. Zero events fall within the signal 
region denned as ±20 MeV from the central value [14| . We set an upper limit of 2.30 events 
in the absence of background at the 90% confidence level (C.L.) for iV@g and 
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FIG. 2: The missing mass distributions of data (left) and Monte Carlo simulation (right) for 
if)(2S) -» K° s pK-n. 



where 0.686 is the decay ratio of Kg to 7r + /T _ 15J; (2.85 ± 0.08)% is the detection efficiency 
and the uncertainty is statistical error of the Monte Carlo sample, and 14.0 x 10 6 is the total 
number of BESII if)(2S) events jlrjj ]. 

Another possibility is that the if>(2S) decays to only one 6 or 6 state. To determine the 
number of 0(1540) events from single 9 or 9 production, we count the number of events 
within regions of 1.52 - 1.56 GeV/c 2 , shown by the arrows in Fig. 3 and set upper limits on 
the branching fractions at the 90% C.L. 
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K° s pK-n) < 0.60 x 10~ 5 

K° s pK + n) < 0.70 x 10~ 5 . 

Backgrounds are not subtracted in the calculation of the upper limits. The numbers used 
to determine the upper limits are summarized in Table I. 

For the decays of J/ip — > K° s pK~n and KgpK + n, we just use the same criteria and 
analysis method as those used for the ip(2S) data to study possible 6(1540) production. 
The scatter plot of K~n {K° s p) versus K^p (K + n) is shown in Fig. 5. The individual 
mass distributions of J /if) data are shown in Fig. 6, which contains 25 KgpK~n events 
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FIG. 3: Mass distributions of i/j(2S)—>KgpK n and KgpK + n modes. 
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FIG. 4: Scatter plot of K n (Kgp) versus K^p (K + n) for ip(2S) — ► K^pK n + KgpK + n modes. 



Decay mode 


N obs 


Efficiency 


Upper limit 


Y>(2S) -» 9(9 -> K° s pK~n 








+K G s pK+n 





(2.85 ±0.08)% 


0.88 x 10~ 5 


i/j{2S) -» Orii -» KgpK~n 


1 


(4.07 ±0.09)% 


1.0 x 10~ 5 


ij){2S) -» GK+n -» K s pK + n 


4 


(3.17 ±0.08)% 


2.6 x 10~ 5 


ip(2S) -» Jfgp0 -» K%pK~n 





(3.99 ±0.09)% 


0.60 x 10~ 5 


ip{2S) -» i^p9 -» K® s pK + n 





(3.42 ±0.08)% 


0.70 x 10~ 5 



TABLE I: Summary of numbers used in the determination of upper limits for the ip(2S) data. 



and 21 K s pK + n events. There is no 6(1540) signal, and we determine upper limits on the 
branching fractions at the 90% C.L. 
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FIG. 5: Scatter plot of if~n (Kgp) versus Kgp (K + n) for J/ijj — > K s pK~n + KgpK + n modes. 



#(J/V> -> 00 -> K%pK~n + K° s pK + n) < 1.1 x 10~ 5 
i3(J/^ -> eK"n -> K s pK-n) < 2.1 x 10~ 5 
#(J/V> -> 0i^ + n -> K° s pK + n) < 5.6 x 10~ 5 
-> Xgp0 -> K s pK-n) < 1.1 x 10~ 5 
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FIG. 6: Individual mass distributions for the J/ip — ► K^pK n and K^pK + n modes. 
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These results along with the numbers used to determine them are summarized in Table II. 



V. SUMMARY 



In this work, we studied if)(2S) and J/ip hadronic decays to KgpK n and KgpK + n 
using 14M ip(2S) and 58M J/ip events. No 0(1540) signal is observed. We set upper 
limits for B(ip(2S) -> 00 -> K% V K~n + K° s pK + n) < 0.84 x lO" 5 , B(J/^ -> 00 -> 
KgpK~n + K° s pK + n) < 1.1 x 10~ 5 at 90% confidence level (C.L.). For the case of single 
0(1540) production, the upper limits determined are also on the order of 10~ 5 for both 
ip(2S) and Jjijj data. 



Decay mode 



N u Efficiency Upper limit 



J/ip -» 96 -> K%pK'n 







1 


(0.88 ± 0.04)% 


1.1 x 10" 


-5 


J/lp- 


-> eif-n -» K^pK-fi 


4 


(0.96 ± 0.04)% 


2.1 x 10" 


-5 


J/lp- 


-> 9if + ra -» K Q s pK + n 


8 


(0.59 ± 0.03)% 


5.6 x 10" 


-5 


J/lp- 


-» Kgj>9 -> K%pK-fi 


2 


(1.19 ±0.05)% 


1.1 x 10" 


-5 


J/v- 


-» K° s p<d -> K^pK+n 


2 


(0.86 ± 0.04)% 


1.6 x 10" 


-5 



TABLE II: Summary of numbers used in the determination of upper limits for the J/ip data. 
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